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We report a measurement of the bottom-strange meson mixing phase /3 S using the time evolu- 
tion of B® — > J/t/>(-> fj, + n~ )<£(->• K + K~) decays in which the quark-flavor content of the bottom- 
strange meson is identified at production. This measurement uses the full data set of proton- 
antiproton collisions at yfs = 1.96 TeV collected by the Collider Detector experiment at the 
Fermilab Tevatron, corresponding to 9.6 fb _1 of integrated luminosity. We report confidence re- 
gions in the two-dimensional space of j3 s and the B® decay- width difference AF S , and measure 
ft s e [— 7r/2, -1.51] U [-0.06,0.30] U [1.26,7r/2] at the 68% confidence level, in agreement with the 
standard model expectation. Assuming the standard model value of /3 S , we also determine Ar s = 
0.068±0.026(stat)±0.009(syst) ps" 1 and the mean B° lifetime, r s = 1.528±0.019(stat)±0.009(syst) 
ps, which are consistent and competitive with determinations by other experiments. 
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The noninvariance of the physics laws under the simul- 
taneous transformations of parity and charge conjugation 
(CP violation) is accommodated in the standard model 
(SM) through the presence of a single irreducible com- 
plex phase in the weak-interaction couplings of quarks. 
A broad class of generic extensions of the SM is expected 
to naturally introduce additional sources of CP violation 
that should be observable, making CP-violation stud- 
ies promising to search for experimental indications of 
new particles or interactions. Thus far, CP violation 
has been established in transitions of strange and bot- 
tom hadrons, with effects consistent with the SM in- 
terpretation [l}[3]. Much less information is available 
for bottom-strange mesons, B° s . Studies of B®-B s fla- 
vor oscillations are unique in that they probe the quark- 
mixing (Cabibbo-Kobayashi-Maskawa, CKM) matrix el- 
ement Vts, which directly enters the mixing amplitude. 
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Large non-SM enhancements of the mixing amplitude are 
excluded by the precise determination of the oscillation 
frequency in 2006 @). However, non-SM particles or cou- 
plings involved in the mixing may also increase the size 
of the observed CP violation by enhancing the mixing 
phase /3 S = arg[— (V ts V t * b ) / (V cs V* b )] with respect to the 
value expected from the CKM hierarchy /3f M sa 0.02 0, 
henceforth referred to as 'SM expectation'. A non-SM 
enhancement of (3 S would also decrease the size of the 
decay-width difference between the light and heavy mass 
eigenstates of the meson, AT S = — T# . The val- 
ues of the mixing phase and width difference are loosely 
constrained, and currently the subject of intense exper- 
imental activity. The analysis of the time evolution of 
B® — > J/00 decays provides the most effective determi- 
nation of f3 s and Ar s [BJ. Assuming negligible contribu- 
tions from sub- leading decay amplitudes [7], the under- 
lying b— > ccs quark transition is dominated by a single 
real amplitude, making j3 s the sole CP-violating phase 
observable, through the interference between the ampli- 
tudes of decays occurring with and without oscillations. 

The first determinations of (3 S , by the CDF and DO 
experiments, suggested a mild deviation from the SM ex- 
pectation [5] . The interest in this measurement increased 
further recently, because of the 3.9ct departure from the 
SM expectation of the dimuon asymmetry observed by 
DO in semileptonic decays of B^ mesons [9], which is 
tightly correlated with /3 S , if generated in the B° sec- 
tor [5] . While updated measurements in £?° — > Jfip<f> de- 
cays [T0HT3] showed increased consistency with the SM, 
more precise experimental information is needed for a 
conclusive interpretation. 

In this Letter we report a measurement of /3 S ; Ar s ; 
the mean lifetime of heavy and light B® mass eigenstates, 
t s = 2/ (V h + Tl); and the angular momentum compo- 
sition of the signal sample using the final data set col- 
lected by the CDF experiment at the Tevatron proton- 
antiproton collider, corresponding to an integrated lumi- 
nosity of 9.6 fb _1 . The analysis closely follows a pre- 
vious measurement in a subset of the present data |10j . 
and introduces an improved determination of the sam- 
ple composition based on a new study of the K + K~ and 
J/ipK + K~ mass distributions. 

The CDF II detector is a magnetic spectrometer sur- 
rounded by electromagnetic and hadronic calorimeters 
and muon detectors that has cylindrical geometry with 
forward-backward symmetry. Charged particle trajecto- 
ries (tracks) are reconstructed using single- and double- 
sided silicon microstrip sensors arranged in seven cylin- 
drical layers [H] and an open cell drift chamber with 96 
layers of sense wires [T5], all immersed in a 1.4 T axial 
magnetic field. The resolution on the momentum compo- 
nent transverse to the beam, px, is a PT /p T « 0.07% (j>t 
in GeV/c), corresponding to a mass resolution of our B® 
signal of about 9 MeV/c 2 . Muons with px > 1.5 GeV/c 



are detected in multiwire drift chambers |16j . A time- 
of-flight detector identifies charged particles with pr < 2 
GeV/c [T7J, complemented by the ionization-energy-loss 
measurement in the drift chamber at higher transverse 
momenta. The combined identification performance cor- 
responds to a separation between charged kaons and pi- 
ons of approximately two Gaussian standard deviations, 
nearly constant in the relevant momentum range. Events 
enriched in J/0— s-/i + /i~ decays are recorded using a low- 
Pt dimuon online selection (trigger) that requires two 
oppositely-charged particles reconstructed in the drift 
chamber matched to muon chamber track segments, with 
a dimuon mass between 2.7 and 4.0 GeV/c 2 . 

In the analysis, two tracks matched to muon pairs are 
required to be consistent with a J/ip— [i~ decay, with 
dimuon mass 3.04 < m MM < 3.14 GeV/c 2 . These are 
combined with another pair of tracks consistent with a 
(f)^K + K- decay, 1.009 < m KK < 1.028 GeV/c 2 , in 
a kinematic fit to a common vertex. A dimuon mass 
constraint to the known J/0 mass [T] improves the B® 
mass resolution. An artificial neural network (NN) classi- 
fier |10j combines multiple discriminating variables into a 
single quantity that statistically separates the signal from 
the dominant background from combinations of real J/0 
decays with random track pairs and a minor component 
of random four-track combinations (both collectively re- 
ferred to as combinatorics) . The NN is trained with sim- 
ulated events for the signal and data from sidebands in 
B° mass, [5.29,5.31] U [5.42,5.45] GeV/c 2 , for the back- 
ground. In decreasing order of discriminating power, the 
input variables to the NN include kinematic quantities, 
muon and hadron particle identification information, and 
vertex fit quality parameters. 

Figure [l] shows the J/iJjK+K- mass distribution from 
the final sample of candidates that pass an NN threshold 
chosen as to maximize the sensitivity to the measurement 
of /3 S [10] . The distribution shows a signal of approxi- 
mately 11 000 decays, above a fairly constant background 
dominated by the prompt combinatorial component, and 
smaller contributions from mis-reconstructed B decays. 

We determine the quantities of interest using a fit to 
the time evolution of bottom-strange mesons. The dif- 
ferences in time evolution of states initially produced as 
a Bg or B s meson are included in the fit as well as the 
differences between decays that result in a CP-odd or 
CP-even combination of the J/00 angular momenta. The 
proper decay time of a B® candidate is a fit observable 
calculated as t — ML xy /pT, where L xy is the distance 
from the primary vertex to the B® decay vertex, pro- 
jected onto the B® momentum in the plane transverse to 
the beam, pr; and M is the known mass of the B® me- 
son pp. The proper decay-time uncertainty, at, is cal- 
culated from the measurement uncertainties in L xy . Be- 
cause the B® meson has spin zero and J/0 and have 
spin one, the B°— > J/ipcf) decay involves three possible an- 
gular momentum states of the J/00 system. These are 
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FIG. 1: (Color online) Distribution of J/ipK + K~ mass with 
fit projection overlaid. 

combined into three polarization amplitudes, longitudi- 
nal polarization (Aq), an d transverse polarization with 
the linear polarization vectors of the vector mesons par- 
allel (Ay) or perpendicular (A±) to each other. The first 
two states are CP even, while the last state is CP odd. 
A CP-odd state can also be produced by a nonresonant 
K + K~ pair or can originate from the decay of the spin- 
/o(980) meson, which results in another independent 
decay amplitude, the S- wave A$. 

To enhance the sensitivity to j3 s , the time-evolution 
of the four decay amplitudes along with six interference 
terms is fitted simultaneously by exploiting differences 
in the distribution of the kaons' and muons' decay an- 
gles. The angles are parametrized in the transversity 
basis, p = (cos0,$,cos\l/) [18], which allows a conve- 
nient separation of the CP-even and CP-odd terms in 
the likelihood. Reference QI5] details the expression for 
the decay rate differential in the decay time and angles. 
The rate is a function of the physics parameters of in- 
terest, /3 S , Ar s , t 8 , and the decay amplitudes with their 
CP-conserving phases. For these we choose A to be real 
and define the CP-conserving phases as Sii = arg (A||/A )> 
S± = arg (A±/A ) and 5s = arg (As/A ). The decay rate 
is also a function of the P|? mixing frequency, which is a 
fit parameter constrained to the experimental value mea- 
sured by CDF, Am s = 17.77 ± 0.12 ps" 1 @]. 

The flavor of the meson at the time of production is 
inferred by two independent classes of flavor tagging al- 
gorithms [TU], which exploit specific features of the in- 
coherent production of bb quarks-pairs in pp collisions. 
Using flavor conservation of the strong interaction, the 
opposite-side flavor tag (OST) infers the signal produc- 
tion flavor from the decay products of the b hadron pro- 
duced by the other b quark in the event by using the 



charge of muons or electrons from semilcptonic B decays 
or the net charge of the opposite-side jet. The same-side 
kaon tag (SSKT) deduces the signal production flavor by 
exploiting charge-flavor correlations of the neighboring 
kaons produced during its fragmentation. The fraction 
of candidates tagged by a combination of OST algorithms 
totals £ost = (92.8 ±0.1)%. The probability of wrongly- 
tagging the meson, wost, is determined per event and 
calibrated using 82 000 B ± J/fy>(-» p + p^)K ± decays 
fully reconstructed in the same sample as the signal [20] . 
Because the P does not oscillate, the OST tag is com- 
pared with the actual flavor, known from the charge of 
the meson. A single scale factor that matches the 
predicted mistag probability to the one observed in data 
is then determined to be 1.085 ±0.035. The observed av- 
eraged dilution, Post = 1 — 2iy ST> equals (12.3 ±0.6)% 
including the scale factor, resulting in a tagging power of 
£ost£>ost = ( L39 ± 0.05)%. The SSKT algorithms tag 
a smaller fraction of candidates, esskt = (52.2 ± 0.7)%, 
with better precision. In the — ¥ J/4>4> sample an aver- 
age dilution of Psskt = (25.9 ±5.4)% is achieved includ- 
ing a 0.94 ± 0.20 scale factor obtained by measuring the 
B® oscillation amplitude in approximately 11 000 (1 850) 
B®— >D~Tr + (Tr + ir~) decays reconstructed in the data cor- 
responding to the first 5.2 fb" 1 [ID]. The resulting SSKT 
tagging power is Ssskt^sskt = (3.5 ± 1.4)%. Higher in- 
stantaneous luminosity conditions in later data resulted 
in a reduced trigger efficiency for hadronic B® decays. 
Hence, the additional sample of P° — > DJtt + (tt + tt^) de- 
cays is too limited for a significant test of the SSKT per- 
formance. Because the SSKT calibration is known for 
early data only, we conservatively restrict its use to the 
events collected in that period. Simulation shows that 
this results in a degradation in (3 a resolution not exceed- 
ing 15%. 

The unbinned maximum likelihood fit uses 9 observ- 
ables from each event to determine 32 parameters in- 
cluding j3 s and Ar, other physics parameters such as P° 
lifetime, amplitudes and phases, and several other quan- 
tities, called nuisance parameters, such as tagging dilu- 
tion scale factors. The fit uses the information of the 
reconstructed B® candidate mass and its uncertainty, m 
and a m ; the B® candidate proper decay time and its un- 
certainty, t and at] the three transversity angles, p; and 
tag information, T> and £; where T> is the event-specific 
dilution given by the mistag probability, and £ is the tag 
decision. Both tagged and untagged events are used in 
the fit. The single-event likelihood is described in terms 
of signal, P s , and background, P b , probability density 
functions (density henceforth) as 

C cx f s P s {m\a m )P s {t,p^\V,(7 t )P s (a t )P s (V) 

+ (1 - f s )P b (m)P b (t\<j t )P b (p)P b (a t )P b (V), (1) 

where f s is the fraction of signal events. The signal mass 
density, P s (m\a m ), is parametrized as a single Gaus- 
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sian with a width determined independently for each 
candidate. The background mass density, p,(m), is 
parametrized as a straight line. The time and angular 
dependence of the signal, P s (t, p,^,\T>,a t ), for a single 
flavor tag are written in terms of two densities, P for B® 
and P for P°, as 

(^^P{t,M) + l ~^P{t 1 p\cj t )^e{p) 1 (2) 

which is extended to the case of OST and SSKT inde- 
pendent flavor tags. Acceptance effects on the transver- 
sity angle distributions are modeled with an empirical 
three-dimensional joint probability density function ex- 
tracted from simulation, e(p), whose largest excursions 
do not exceed 15%. The time and angular distributions 
for flavor-tagged B® (B®) decays, P (P), are given by 
the normalized decay rate as functions of decay time 
and transversity angles of Ref. |19j . assuming no direct 
CP violation in the decay. Building on previous mea- 
surements |2T], we model the decay-time density for the 
background, Pb(t\a t ), with a ^-function at t = 0, one pos- 
itive, and two negative exponential functions. All time- 
dependent terms are convolved with a proper time reso- 
lution function, modeled as a sum of two Gaussians with 
common mean and independent widths determined by 
the fit. The resulting decay-time resolution is equivalent 
to that of a Gaussian distribution with 90 fs standard 
deviation. The background angular probability density 
is determined from P° mass sideband events to factorize 
as P b (p) = P b (cos9)P 6 ($)P b (cos*). The distributions 
of the decay-time uncertainty and the event-specific di- 
lution differ for signal and background events, thus their 
densities are explicitly included in the likelihood. The 
probability density functions of the decay-time uncertain- 
ties, P s (at) and p,((Tt), are described with an empirical 
model from an independent fit to the data. The signal 
density, P S (T>), is determined from binned background- 
subtracted signal distributions, while the background 
density, Pf,(£>), is modeled from candidates in the signal 
sidebands. Potential sources of systematic uncertainties, 
associated with imprecisely known calibration factors of 
tagging dilutions, are taken into account by floating these 
factors in the fit within Gaussian constraints. 



The likelihood function shows two equivalent global 
maxima, corresponding to the solutions with positive and 
negative value of AT S , and additional local maxima gen- 
erated by approximate symmetries (T5]. Multiple solu- 
tions make the estimation of parameters and their un- 
certainties challenging with limited sample size. If j3 s is 
fixed to its SM value, the fit shows unbiased estimates 
and Gaussian uncertainties for AT s , t s , polarization am- 



plitudes, and the phase 5±, yielding 

t s = 1.528 ± 0.019(stat) ± 0.009(syst) ps, 
Ar s = 0.068 ± 0.026(stat) ± 0.009(syst) ps" 1 , 
|A | 2 = 0.512 ± 0.012(stat) ± 0.018(syst), 
\A\\ | 2 = 0.229 ± 0.010(stat) ± 0.014(syst), 

S ± = 2.79 ± 0.53(stat) ± 0.15(syst). 

The correlation between r s and Ar s is 0.52. We do 
not report a measurement of $|| . The fit determines 
$H ~ 7r, but the estimate is biased and its uncertainty 
is non-Gaussian because the likelihood symmetry under 
the $|| — > 2ir — 8\\ transformation |19j results in multiple 
maxima in the vicinity of $|| = n. Systematic uncer- 
tainties include mismodeling of the signal mass model, 
decay-time resolution, acceptance description, and an- 
gular distribution of the background; an 8% contamina- 
tion by B ^ J/ipK* (892)° and P ^ J/tpK+TT- decays 
misreconstructed as P° — > J/ip(f> decays; and silicon de- 
tector misalignment. For each source, uncertainties are 
determined by comparing the fit results from simulated 
samples in which the systematic effect is introduced in 
the model and samples simulated according to the de- 
fault model. The uncertainty on the Ar s measurement 
is dominated by the mismodeling of the background de- 
cay time. The largest contribution to the uncertainty 
on t s is the effect of silicon detector misalignment. The 
angular acceptance model dominates the systematic un- 
certainties on the amplitudes. 

If Ps is free to float in the fit, tests in statistical trials 
show that the maximum likelihood estimate is biased for 
the parameters of interest, and the biases depend on the 
true values of the parameters. Hence, we determine con- 
fidence regions in the j3 s and (j3 s , AT S ) spaces by using a 
profile-likelihood ratio statistic as a \ 2 variable and con- 
sidering all other likelihood variables as nuisance param- 
eters. The profile-likelihood ratio distributions observed 
in simulations deviate from the expected \ 2 distribution, 
yielding confidence regions that contain the true values of 
the parameters with lower probability than the nominal 
confidence level. In addition, the profile-likelihood ratio 
distribution depends on the true values of the unknown 
nuisance parameters. We use a large number of statisti- 
cal trials to derive the profile- likelihood ratio distribution 
of our data. The effect of nuisance parameters is ac- 
counted for by randomly sampling their 30-dimensional 
space within ha of their estimates in data and using the 
most conservative of the resulting profile-likelihood ra- 
tio distributions to derive the final confidence regions. 
This procedure ensures that the confidence regions have 
nominal statistical coverage whatever the configuration 
of nuisance parameters values and increases the size of the 
P s confidence interval by about 40%. We determine the 
confidence level for 32 x 48 evenly spaced points in P s <E 
[— 7r/2,7r/2] and Ar s e [-0.3,0.3] ps" 1 and smoothly 
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FIG. 2: (Color online) Confidence regions at the 68% (solid) 
and 95% C.L. (dashed) in the (/3 S , AF S ) plane (main panel). 
The standard model prediction is shown as a circle with er- 
ror bars [55]. The inset shows the coverage-corrected profile- 
likelihood ratio as a function of j3 a , in which AF S is treated 
as all other nuisance parameters. 



interpolate between them to obtain a continuous region 
(Fig. [2]) . Assuming the standard model values for j3 s and 
Ar s , the probability to observe a profile- likelihood ra- 
tio equal to or higher than observed in data is 54%. By 
treating Ar s as a nuisance parameter, we also obtain 
(5 S € [-71-/2, -F51] U [-0.06, 0.30] U [F26, tt/2] at the 68% 
C.L., and /3 S e [-tt/2, -1.36]U[-0.21, 0.53]U[1.04, tt/2] at 
the 95% C.L. The fraction of S-wave in the K + K~ mass 
range 1.009-1.028 GeV/c 2 is determined from the angular 
information to be consistent with zero with 0(2%) un- 
certainty, which is in agreement with our previous deter- 
mination [TO] and the LHCb and ATLAS results [T21IF3] . 
and inconsistent with the DO determination llj. An 
auxiliary simultaneous fit of the K + K~ and J/ij)K + K~ 
mass distributions 23J that includes the full resonance 
structure of the J5°— > J/%pK + ir~ decay [H] is performed. 
The K + K~ mass is fit in a range enlarged to 0.988- 
1.2 GeV/c 2 using a relativistic Breit-Wigner distribution 
for the 4> meson, the shape suggested in Ref. [25] for 
the /o(980) meson, and an empiric shape determined 
from data for the combinatorial background. In the 
1.009-1.028 GeV/c 2 mass range, this fit determines a 
(0.8 ± 0.2(stat))% K + K~ S-wave contribution in agree- 
ment with the central fit, and a contamination from mis- 
identified B° decays of (8.0±0.2(stat))%, which is signif- 
icantly larger than the 1-2% values typically derived as- 
suming only P-wave B° decays [TUllllj . If neglected, this 



additional B° component could mimic a larger K + K S- 
wave than present. 

In summary we report the final CDF results on the 
B® mixing phase and decay width difference from the 
time-evolution of flavor-tagged B®—> J/ip(j> decays recon- 
structed in the full Tevatron Run II data set. This analy- 
sis improves and supersedes the previous CDF measure- 
ment obtained in a subset of the present data [TU] ■ Con- 
sidering Ar s as a nuisance parameter, and using the 
recent determination of the sign of Ar s [26) . we find 
-0.06 < P a < 0.30 at the 68% C.L. Assuming a SM value 
for P s , we also report precise measurements of decay- 
width difference, Ar s = 0.068 ±0.026(stat)±0.009(syst) 
ps _1 , and mean B° lifetime, t s = 1.528 ± 0.019(stat) ± 
0.009(syst) ps. All results are consistent with expec- 
tations and with determinations of the same quantities 
from other experiments |11H13] , and significantly improve 
the knowledge of the phenomenology on CP violation in 
B® mixing. 
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